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ABSTRACT 



We investigate an in-situ formation scenario for Earth-mass terrestrial plan- 
ets in short-period, potentially habitable orbits around low-mass stars {M^: < 
O.SMq). We then investigate the feasibility of detecting these Earth-sized plan- 
ets. We find that such objects can feasibly be detected by a ground-based tran- 
sit survey if their formation frequency is high and if correlated noise can be 
controlled to sub-milli-magnitude levels. Our simulations of terrestrial planet 
formation follow the growth of planetary embryos in an annular region span- 
ning 0.036 AU < a < 0.4 AU around a fiducial M7 (O.I2M0) primary Ini- 
tial distributions of planetary er abryos are c alcula ted using the semi-analytic 
evolutionary model outlined by I Chambers! (120061 ). This model specifies how 



planetary embryos grow to the stage where the largest embryo masses lie in 
the lO'^^g < Membryo < 5 X lO"^^ g range (corresponding to the close of the so- 
called oligarchic growth phase). We then model the final phases of terrestrial 
planet assembly by allowing the embryos to interact with one another via a 
full N-body integration using the Mercury code. The final planetary system 
configurations produced in the simulations generally consist of 3-5 planets 
with masses of order 0.1 — l.OM® in or near the habitable zone of the star. We 
explore a range of disk masses (O.2M0 to 3.3M^) to illuminate the role disk 
mass plays in our results. With a high occurrence fraction or fortunate align- 
ments, transits by the planet formed in our simulations could be marginally 
detected with modest telescopes of aperture Im or smaller around the near- 
est M-dwarf stars. To obtain a concrete estimate of the detectability of the 
planets arising in our simulations, we present a detailed Monte-Carlo transit 
detection simulation incorporating sky observability, local weather, a target 
list of around 200 nearby M-dwarfs, and a comprehensive photometric noise 
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model. We adopt a baseline 1.5 mmag level of correlated stellar noise sampled 
from the photometry of the planet-bearing red dwarf Gl 436. With this noise 
model we find that detection of IR^ planets around the local M-dwarfs is chal- 
lenging for a ground-based photometric search, but that detection of planets 
of larger radius is a distinct possibility. The detection of Earth-sized planets 
is straightforward, however, with an all-sky survey by a low-cost satellite mis- 
sion. Given a reduced correlated noise level of 0.45 mmag and an intermediate 
planetary ice-mass fraction of planets orbiting a target list drawn from the 
nearest late-type M dwarfs, a ground-based photometric search could detect, 
on average, 0.3 of these planets within two years and another 0.5 over an in- 
definitely extended search. A space-based photometric search (similar to the 
TESS mission) should discover ~ 17 of these Earth-sized planets during it's 
two year survey, with an assumed occurrence fraction of 28%. 



Keywords: PLANETARY FORMATION ; EARTH ; TERRESTRIAL PLAN- 
ETS ; EXTRASOLAR PLANETS ; TRANSIT 
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1 Introduction 



To date, planets have been detected in every region of mass and orbital pa- 
rameter space to which observational searches are sensitive. In recent years, 
the minimum detected planetary mass has been decreasing at a very rapid 
pace. If one projects this trend forward, one sees that detection of the first 
near Earth-mass extrasolar planets orbiting main sequence stars should oc- 
cur sometime around 2011. Thus, while it has always been of great interest 
to search for terrestrial mass planets, the time when success might be ob- 
tained now seems close at hand. Indeed, both radial velocity observations and 
photometric transit surveys have the capability of achieving t he detection o 
i -mass planets orbiting nearby low -mass stars, (see e.g. 



Eart 



2008 



Guedes et al. 



Nutzman and Charbonneaul . 



20081) 



In the standard paradigm of Solar system formation, the terrestrial planets 
emerge during the later stages of proto-planet ary di sk evolut i on th rough the 



coagulation of smaller bodies (ISafronov 



19691 ) (see 



Lissauerl . 



1993 



tailed review). The initial formation process proceed s 



1993 



Weidenschilling et al. 



1997 



Kokubo and Ida. 



rapid 



1998, 



lie. 



, for a de- 



Wetherill and Stewart 



2000) with a few of 



the small bodies accreting more matter than their companions, thereby in- 
creasing th eir gravitational cross-section s and experiencing even more rapid 



growth (e.g. 



Wetherill and Stewart 



19891 ) This "runaway growth" phase rapidly 



produces planetary embryos from the sea of planetesimals. As the embryos 
grow, they perturb the orbits of local planetesimals, thus reducing gravita- 
tional focusing and slowing the continued growth of the planetary embryos. 
Nevertheless, the largest planetary embryos still accrete material more rapidly 
than any remaining low-mass planetesimals - leading to the so-called "oli- 
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garchic" growth phase. Eventually, an emergent distribution of hundreds of 
planetary embryos interact and collide stochastically, until all the planetary 
embryos have either been incorporated into the forming planets or ejected 
from the system. 

N-body simulations of the final accretion phases of terrestrial planet assembly 



are able to recreate the 
well (e.g. lO'Brien et al. 



jroad -brush aspects of the inner solar system fairly 



20061 ). In particular, the orbital properties of Earth 
and Venus analogs are often reproduced, whereas Mars-like and Mercury-like 
objects are harder to account for. Ultimately, however, one would like to test 
the predictive outcome of the current paradigmatic model on systems other 
than our own. It has been shown, for example, that initial distributions of 
lunar-sized embryos similar to the distributions that can generate the solar 
system ter restrial planets predict Earth- mass worlds o rbiting both Alpha Cen 



A and B (j Quint ana et al. 



2002 



Guedes et al. 



20081). The accretio n of such 



2008 



20091 ) 



embryos, however, may be problematic (iThebault et al.l . 
Accretion simulations can also be extended to investigate the formation of 
planets around the low- mass stars that const i tute t he bulk of the local stel- 
lar population. For example. 



Raymond et al. 



(120071 ) present a suite of accre- 



tion calculations which predict that Earth-sized terrestrial planets should be 
rare around M-type stars with < O.SM©. Their models start with initial 



planetesimal surface densities that scale generally downwarc 



mu m mass solar nebula-b ased models studied by 



the 



Raymond et al 



Chambers 



Tom the mini- 



(1200 ih . Indeed, 



(120071 ) result can be encapsulated by stating that low-mass 



planetesimal disks produce low-mass planets. Their work makes it clear that 
if the protoplanetary disks associated with the lowest mass stars are scaled 
down versions of the Sun's protoplanetary disks, then habitable planets orbit- 
ing red-dwarfs will be very rare. The Raymond et al. result emphasizes the 
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importance of observations probing inner disk masses as the available mass in 
these regions is pivotal to the formation of habitable planets in-situ. 
Planetary formation around M-dwarfs is now the focus of considerable ob- 
servational and theoretical work. One conclusion that seems robust is that 



Jovian-mass gas giant planets aroun d M-dwarfs sho u. 
masses scale with stellar mass (e.g. iLaughlin et al 



d be rare if outer disk 



200J). This prediction 



is a natural consequence of the core ac c retion theory for giant planets (see, 



e.g. 



Pollack et all. 119961: IH 



observations ( Butler et al. 



jickyi et al. 



2004 



2005) 



EndletaL 



, and now has supporte d from 



2006 



Lovis et al. 



2006|). Fur- 



thermore, ice-giant planets akin to Uranus and Neptune should be common 
around M dwarfs. This result is lent plausibility by the radial velocity obser- 
vations of 



Mayor et al 



(120081 ) who estimate the fraction of FGK s tars with 



ice-gi ant planets to be at least 30%, and by microlensing results (e.g. 



Bennett 



20091 ) which give direct evidence of such planets. Since accretion models are 
so dependent on the available mass, the seeming plethora of planets is hard 
to understand if the only available initial condition is a scaled version of the 
MMSN. Thus, because M stars allow for ready detection of companions down 
into the terrestrial mass range, we think it is of interest to investigate the 
observable consequences of a broad range of possible formation scenarios. 
The plan for this paper is as follows. In §2, we describe the details of our terres- 
trial planet formation model for low-mass stars. This model is designed from 
the outset to favor the production of habitable planets orbiting red dwarfs. 
In §3, we describe the specific simulation outcomes for our model, and aggre- 
gate their statistical properties. In §4, we construct end-to-end Monte-Carlo 
simulations that realistically assess the near-term prospects for finding these 
planets using inexpensive photometric surveys that target nearby stars. In §5, 



we conclude with our finding that a ground-based survey with small dedicated 
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telescopes has difficulty providing significant constraints on our scenario, but 
that an inexpensive microsatellite survey can provide a significant test. 



2 Formation Simulation Method 



We have carried out 40 N-body accretion simulations of terrestrial planet for- 
mation around low mass stars. The simulations use initial conditions that are 
generated with a semi-analytic model that is very similar to that described in 



Chambers! (120061 ). The semi-analytic evolutionary model is intended to encap- 
sulate the growth of bodies up to and including the oligarchic growth phase, 
and includes (1) a disk of equal-mass planetesimals with their equilibrium 
eccentricities and inclinations, (2) a gas disk which acts to damp these ec- 
centricities and inclinations, and (3) seed planetary embryos which accrete 
material throughout the simulation and dynamically heat the planetesimal 
distribution. As the planetary embryos grow, they accrete and retain atmo- 
spheres, which increase their collisional cross-sections, and further accelerate 
their growth. A semi-analytic evolution is deemed complete when the distribu- 
tion of bodies in the model reaches 50 planetary embryos (each separated by 
10 mutual Hill Radii). At this point the mass in planetary embryos and plan- 
etesimals is approximately equal (~ 1.6M^ in each population), and we take 
the orbital parameters of the embryos and the surface density profile of the re- 
maining planetesimals and fragments as the initial condition and transfer the 
system to a 3 dimensional N-body simulation. Our N-body simulations begin 
with the 50 planetary embryos and include 50 additional particles as represen- 
tatives of the remnant planetesimal population. Because we begin with equal 
numbers, the masses of the 'super-planetesimals' and the planetary-embryos 
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are approximately equal (although as the simulations proceed, the embryos 
rapidly outpace the planetesimals) . 

The bodies orbit a late type M-dwarf star (M7, = 0.12 Mq) and are 
endowed with random orbital eccentricities in the range [0.00 < Cj < 0.01], 
and random orbital inclinations in the range [—0.5° < ii < 0.5°]. The initial 
longitudes of periastron, zui, longitudes of the ascending nodes, Qi, and mean 
anomalies. Mi, of the bodies are chosen randomly. 

Due to the uncertainties in the surface density profile and to facilitate com- 
parison with previous work, our semi-analytic model draws on an initial solid 
surface density profile chosen to resemble the Minimum Mass Solar Nebula 



(IWeidenschilling 



19771 ) in that it falls as S = an • (a/l AU) '^^^ unti l it reaches 



2004 ) where we 



the 150 K ice-sublimation temperature (iPodolak and Zuckerl . 
truncate the outer edge of the disk. This choice of surface density profile plays 
a significant role in our investigation - given a shallower profile, formation of 
terrestrial planet s in M -dwarf habitable zones in-situ becomes more difficult 



([Raymond et al 



20071 ). and under these conditions, the smaller planetary 
population would be undetectable given current g round-based photom etric 



precision . In o ur model, spurred on the results of 



Bennett 



Mayor et al 



fl2008h and 



( l2009l ) we imagine that an ice-rich (and up to Neptune-mass) object 
has already formed in the region exterior to the ice-line due to the larger initial 
solid surface densities there (in analogy with the formation of Jupiter prior 
to the formation of the Earth). Given our fiducial luminosity = 0.02Lq, 
the 150 K fiducial ice-line lies at risox = 2.7 AUy^L^/Lo = 0.4 AU. We note 
that this is only an approximate boundary. In order to account for the de- 
structive nature of collisions in the inner regions (where coUisional speeds are 
faster) we truncate the inner edge of the disk by taking Chambers' (2001) 
inner disk Keplerian cutoff velocity: v = 54km/s at = 0.3 AU. Around our 
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low-mass primary (M* = O.I2M0) this velocity sets the inner edge of the disk 

to Tinner = 0.036 AU. 



The value of an 



Raymond et al. 



'or th e minimum mass solar nebula is ao^MMSN = 7g/cm . 
(120071 ) show unambiguously that when the MMSN is extrap- 
olated to lower disk masses, the resulting terrestrial planets rarely exceed the 
mass of Mars. We thus ask: is there any evidence beyond the MMSN to which 
we can look that might suggest higher surface density normalizations for the 
disks orbiting the lowest mass stars? 



Sub-m illimeter flux measurements of low- mass primary systems (e.g. lAndrews and Williams 



20071 ) suggest a possible linear correlation between stellar-mass and disk-mass 
for low-mass stars, albeit with a large amount of scatter. To quantify this 
effect, we calculate occurrence fractions of Earth- mass or larger circumstel- 
lar disks given a mean solar type disk-mass of O.OOSM© ± O.bdex and scaling 
down to our stars' disk masses by assuming that the disk-mass scales with the 
stellar-mass to the 0th, 1st or 2nd power. For these three scaling possibilities 
(0, 1, 2) we find occurrence fractions of (28%, 0.7%, 0.001%), respectively. 
Further, there is a scatter of around an order of magnitude on either side of 
the MMSN value of Mdi.sk ~ O.OIM0 for a given stellar-mass primary. This 
serves to restrict the likely surface-density normalization to between 0.1 and 
10 times the MMSN value of ctq^mmsn = 7g/cm^. 

Searching further afield, we can draw clues from other satellite systems or- 
biting primaries other than the Sun. First, consider the Jovian satellite lo: 
Jupiter's mass is O.OOIM0 and lo, with P = 1.8 d, has M = 8.93 x lO^^g. 
If we spread lo's mass out following the r~^/^ surface density profile, we find 
S ~ ll,000g/cm^ at lo's orbital radius. Next, consider the GJ 876 system. 
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20051 ). Planet 'd' has a P = 2 day orbit 



with = O.32M0 fiRivera et all , 
(a = 0.02 AU) and M = 7.5M®. If one assumes that planet d formed from 
material extending to r ~ 0.075 AU and if this material is distributed as 
S oc a~^/^, then we again find S ~ 1.1 x 10"^g/cm^ at the P = 2 d orbital 
radius. 

Going out on a limb, one could argue that this similarity between lo and 
GJ 876 d implies that low-mass primaries can have S ~ 1.1 x lO^g/cm^ 
at the 2-day orbital radius. We adopt this surface density normalization in 
our simulations. Certainly, such an assumption is, at this point, mere conjec- 
ture, but by the same token, a large downward extrapolation from the MMSN 
seems equally uncertain. In our model, if we extrapolate to 1 AU, one finds 
o"o = 21g/cm^. This value, while high, does fall within the large range of ac- 
ceptable surface densitie s inferred from the extant sub-millimeter data. 



Canup and Ward! (120061 ) have noted the striking similarity in mass ratios be- 
tween the Jovian planets (Jupiter, Saturn, Uranus) and their regular satellite 
systems: rris/Mp ~ 2 x 10~^. They suggest that this points to a common for- 
mation scenario within gas-starved circumplanetary disks. In our red dwarf 
context, the presence of an exterior ice-giant planet may play a similar role. 
A Neptune- mass planet orbiting a O.IM0 primary at these radii could me- 
ter the inward fiow of gas through the gap that it has opened in the proto- 
planetary disk. The details of the infiow are dependent on the details of the 
gas fiow around the co-rotational resonances and so is largely undetermined 
for the present situation. We hypothesize that the slow metering of gas into 
the inner region allows planetesim als to build up wit h out s uffering the catas- 



trophic effects of Type I migration 



Canup and WardI (120061 ). This would lead 



to a reduced gas-density in the interior orbital regions which would signifi- 
cantly reduce the effect of Type I migration. In the present work we make 
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the assumption that such a mechanism makes the effects of Type I migration 
neghgible over the planet mass range studied, but we have not done a detailed 
analysis and this scenario would definitely require further investigation should 
observations reveal the presence of a terrestrial-mass planetary population or- 
biting nearby M dwarfs. 

Terrestrial planets that form in a mode analogous to Jovian satellites would be 
expected to be relatively volatile-poor, while volatile-rich planets (which would 
have larger radii at a given mass) w ould be indicative o f a planet that migrated 



inward rather than forming in-situ (IGaidos et al. 



2007 



Raymond et al. 



20081 ). 



The lowest- mass M dwarf target s tars are near 



) y, and so follow-up radial ve- 



locity measurements are possible (iBouchy et al 



20051 ). An Earth mass planet 



in a 30 day orbit around our fiducial O.I2M0 star would elicit a radial velocity 
half-amplitude oi K = 0.85m/s and so is detectable given current technology 



( Lovis et al. 



20081 ). Measurement of the Doppler wobble generates a mass, and 
when combined with the radius from the transit photometry could distinguish 
between the two scenarios. (This situation could, of course, be complicated 
by the presen ce of signi ficant gas atmospheres augmenting the planetary radii 



Adams et al. 



[e.g. 



20081 )). Assuming similar volatile distributions to solar-type 



stars, it seems likely that water content levels o n these planets would be simi 



lar to those predicted by 



Raymond et al. 



(I2OO4J ). However, for the purposes of 



probing bulk composition, transits can only make zeroth-order comparisons. 
It is interesting to speculate as to the water content of these terrestrial planets 
but, as yet, the situation is largely unconstrained. 

Given these considerations, we set the initial solid surface density around our 
M7 star to be S = 21 x (y^) g/cm^, where almost all of this mass is in 
the form of a population of 5 km (10^^ g) planetesimals. The disk scale-height 
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is H = 0.05 iij TjJ ^AU. We adopt constant-mass planetesimals, in keeping 

with iChamberd (120061 ) . This assumption greatly simplifies the semi-analytic 
framework with the drawback of possibly introducing a modest overestima- 
tion of the embryo growth rates. The remainder of the solid surface density 
(0.03% by mass) i s allotted to plan e tary e mbryos in accordance with the con- 



straint, derived by 



Ida and Makind (119931 ) , that the onset of oligarchic growth 



occurs when the embryo and planetesimal masses and surface-densities satisfy 
2Memb^emb = mpT,^,. Wc Set the aver age inter-embryo spacing to be 10 Rhui in 



keeping with 



Kokubo and Idal (119961 ). We also include a gas disk which serves 



to damp the planetesimal's eccentricities and inclinations and causes the plan- 
etesimals to drift inwar ds via gas drag. We implement the effects of gas drag 



using the prescription of 



Adachi et al. 



(119761 ) . The gas disk begins with an ini- 



tial surface density of 200 times the solid surface density and, over the course of 
the simulation, the gas dissipates according to T^gasit) = Egas(O) exp(— t/r^eb), 
where Tmb = 2 Myr. 

As the random velocities of the planetesimals rise, collisions between planetes- 
imals will become destructive, giving rise to many small collision fragments. 
Because these fragments are so small, the gas is very effective at damping out 
the eccentricities and inclinations of these fragments. Thus we treat this pop- 
ulation as having zero eccentricities and inclinations. This greatly strength- 
ens gravitational focusing between the embryos and the fragments, thereby 
increasing the growth rate of embryos substantially over what would occur 
without fragmentation. We also include the effects of the embryo's nascent 
planetary atmospheres. This also tends to increase the embryo's collisional 
cross-sections, further increasing the growth-rates. With all of these factors 
included we find that our embyro growth-rates are a few (~ 3) times faster 
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Chambers! (120061) ( eqn. 6 ), and 



Kokubo and Idal (120021 ) (eqn 



than the equihbrium-regime approximation of 
~ 8 times faster than the estimation given by 
15). 

We run the semi-analytic model until it produces 50 planetary embryos sepa- 
rated by 10 Hill Radii (at which time the embryos have grown to sizes between 
1 X 10^^ and 5 x 10^^ grams). We then replace the remaining planetesimal sur- 
face distribution with 50 "super-planetesimals" (each with nip ~ 2 x 10^^ g). 
The 50 embryos and the 50 super-planetesimals are e volved using the Mercury 



integrator package (jChambers and Migliorini 



19971 ) which employs a hybrid 



algorithm that integrates the bodies with a symplectic map and switches to 
Bulirsh-Stoer integration when close encounters occur. In this new framework, 
the embryos are self-gravitating while the super-planetesimals are not. Each 
N-body simulation is run for 100 Myr to help ensure the long term stability of 
the emergent systems. The hybrid integrator uses a variable stepsize to meet 
a Bulirsch-Stoer step-wise error tolerance of 10^^^. The typical energy error 
at the end of one of our integrations is ~ 4 ■ 10~^. 

In order to address the uncertainty in the disk mass used we have run sim- 
ulations with different disk surface-density normalizations. If the total disk- 

I ^rn 

mass scales as the stellar mass as is implied by lAndrews and Williamsl (120071 ) 



and if the disk surface density normalization in our region scales as the to- 
tal disk-mass, then the mean normalization for the disk surface-density, (Tq ~ 
7g/cm^ X 0.12 = 0.84 g/cm^. The resulting inner disk would only contain 
O.17M0 and thus cannot host detectable terrestrial planets formed in situ. 
In order to quantify this situation, we have performed simulations of plan- 
etary formation and detection with normalizations of ctq = 0.84 g/cm^ and 
(To = 7g/cm^. The resulting detection statistics are low to zero, as expected, 
and are presented for completeness in Table 2. 
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3 Formation Simulation Results 



Our semi-analytic model seeks to simulate the final, rapid, stages of oligarchic 
growth that follow the runaway growth phase. Initial planetary embryo masses 
are Mem6 ~ 3 x 10^^ g. The semi-analytic model allows the planetary embryos 
to achieve presumably more natural mass and spatial distributions, leading to 
a more robust set of initial conditions for the final, stochastic, phase of growth. 
In Fig. 1, the thin and thick lines represent beginning and end states of the 
semi-analytic phase of the simulation. The solid surface density begins al- 
most entirely in the form of planetesimals and after 100 years of growth, 
the inner region has become dominated by the planetary embryos. At this 
point, the planetary embryos are l arge enough that, maintaining the average 



inter-embryo spacing suggested by 



Kokubo and Ida 



(119981 ). the surface den- 



sity profile can support 50 such embryos. Using simple time scale analysis, we 
can compare this oligarchic growth rate with other calculations. We find that 
Tgrow obtained at a = 0.1 AU around our M7 dwarf is 29 times shorter than 
that found in sim ulations designed to mimic the solar system at a = 1 AU 



( 1 Chambers 



20061 ) . Given the approximate nature of timescales, this is consis- 



tent with the acceleration to be expected (see §3.3 for a detailed calculation). 



Once the semi-analytic model has generated initial conditions, we transfer the 
remaining planetary embryos to the N-body simulations. At this transition, 
the gas is suddenly removed. Figures 2 and 3 show the evolution of two of the 
simulations, one from Group A (with nothing exterior to the ice-line), and one 
from Group B (with a Neptune-mass body on a 0.6 AU nearly circular orbit). 
The sizes of each of the symbols are proportional to the radius of the body. 
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Figure 2 shows the evolution of simulation A-01. After 3,000 years, mutual 
interactions have slightly increased the eccentricities of the bodies, and we see 
the rapid buildup into protoplanets of approximately Mars-mass. By 0.1 Myr, 
we see the formation of a few bodies with Membryo — 0.7M®, predominantly 
in the innermost regions of the disk where the surface density is high and 
where the timescales are short. These massive embryos rapidly clear out the 
surrounding regions of lower mass planetary embryos through collisions and 
scattering. Over a timescale of order r ~ 1 Myr, these larger bodies slowly 
scatter or accrete the remaining planetary embryos. Finally, those larger re- 
maining bodies settle into low eccentricity orbits through the ejection or ac- 
cretion of the remaining smaller planetary embryos, producing a system of 
planets with masses between Mars' and Earth's. These remaining planets are 
in low eccentricity orbits separated by, on average, 30 mutual Hill radii. Thus, 
while long term stability is not assured for these s ystems, it is quite likely for 



extremely long timescales (IChambers et al 



19961) 



Figure 3 shows the evolution of simulation B-01. (the Group B simulations 
included a Neptune-mass planet at 0.6 AU) The external perturber's influence 
is noticeable early on through the excitation of eccentricities of several plan- 
etary embryos around the 2:1 resonance at 0.378 AU as well as the increased 
eccentricities of the entire population relative to that seen in the early evo- 
lution of simulation A-01 (Fig. 2). The orbits rapidly become excited due to 
mutual interactions amongst themselves and the perturber, resulting in a dy- 
namically warm system. Within approximately 10,000 years, the system has 
built up a handful of Mars-mass planets through collisions between the plan- 
etary embryos. The growth rate then slows considerably and over the later 
stages of the disk's evolution the larger surviving planets eject most of the 
other remaining low mass bodies with the help of the massive perturber. The 



16 



surviving planets are on low-eccentricity orbits which are separated by, on 
average, 17 mutual Hill radii. The overall situation is similar to that of Group 
A: the planetary embryos have been efficiently ejected or accreted onto the 
remaining planetary bodies which are between Mars' and Earth's mass and 
whose orbits are located from the star's habitable zone outwards. 

An accretion timescale of 10,000 years seems surprisingly short, but it is nev- 
ertheless a natural consequence of the fact that the growth rate scales as 
'Tgrowth oc V ■ . Although our region of interest (a ~ 0.1 AU) is close to the pri- 
mary, the star is also less massive {M^ = O.I2M0) resulting in very similar or- 
bital velocities to the solar system [v ~ l.lf®). For this planet-forming region 
we can compare the local density to the solar system; p/pmmsn = 4.67, giving 
our growth timescale relative to the solar system as: Tgrowth/ — 24. Thus our 



simulations after 100,000 years of evolution should look sim ilar to simu 



ations 



2OO1I). 



done of the solar system at approximately 2.4 million years (iChambersl . 
This is indeed the case. 



Considering planets that exceed the mass of Mars by the end of our simu- 
lations, we see that their mass accretion is 50% complete at t^q = 8 x 10^ 
yr and 90% complete at rgo = 10^ yr. The rapid initial evolution arises from 
the initial flurry of accretion activity in the dense environments. The presence 
of an external perturber accelerates this timescale significantly: the 50% and 
90% stages are reached at median times of t^q = 2x 10* yr and rgo = 0.2 x 10^ 
yr, respectively. 



With 20 complete simulations in both groups A and B, we can consider the 
statistics of our systems. Figure 4 shows the resulting histogram of the plane- 
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tary masses. The two distributions are roughly similar over most of the mass 
range, implying a concordance in the average number and masses of plan- 
ets. At the lowest masses, however, we see the primary effect of an external 
perturber on the surviving populations. The ice-giant is very efficient at eject- 
ing low mass planetary embryos, thus the final planetary mass distribution 
of Group B is made up of a population which is distributed around O.TM^. 
In addition. Group A has a low-mass remnant population that survived with 
few if any collisions and only mild scattering events. We wish to isolate the 
general properties of the two distinct populations of surviving bodies. We thus 
draw a division between the two at 0.05M®. The resulting averaged quantities 
of interest are presented in Table 1. The effect of inclusion of an ice-giant is 
that it serves to raise eccentricities of the low-mass planetary embryos from 
the simulations leading to ejection or accretion onto the star. This decreases 
the average planetary mass slightly, and also limits the growth of the remain- 
ing planets to the inner disk regions. Previous wor k that compared terrestria l 



planet formation with and without giant planets (ILevison and Agnorl . 



20031) 



found that when giant planets caused significant mass-loss from the system 
the resulting terrestrial planets tended to be smaller and closer in to their host 
stars. This is consistent with our findings: when the ice-giant was included we 
have significant mass-loss and find that our resulting planets are smaller and 
closer in to the star. Figure 5 shows the mass found in each bin of semi-major 
axis for Groups A and B. Note that the majority of the massive planets are lo- 
cated inside or very close to the star's so-called habitable zone (0.03 — 0.08 AU) 



(jSelsis et al. 



20071 ). In any case, planetary habitability depends on a variety 
of factors, and in general can only be evaluated once the precise orbit, stellar 
environment and bulk structure of a given candidate have been determined. 
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4 Detection Simulation 



Our simulations create potentially detectable planets. Our formation scenario 
is therefore testable, since a M7 (~ O.145i?0) primary presents a generous 
0.4% transit depth for a terrestrial (~ 1-Re) planet. Furthermore, habitable 
planets orbiting such stars enjoy a relatively high a-priori transit probabilities 
in comparison to an Earth-analog orbiting a solar type star. The a-priori geo- 
metric transit probability for a habitable planet emerging from our simulations 
is of order 

P_., ^ 0.0045 ( ^) ^+-7'-^,^--> . 1%, (1) 
a \ Kq ) 1 — 

where w is the longitude of periastron referenced to the plane of the sky, and 
the other symbols have their usual meanings. 

Given a transit depth and a specified detection strategy, we can compute 
the detectability of our model planets. Figure 6 shows the average detection 
probability per stellar system observed as a function of detection sensitivity 
threshold. For example, if the noise can be controlled to allow the detection of 
planets generating central transit depths deeper than 0.4%, then each planet- 
bearing star yields a P ~ 3.2% chance of a detection. The situation thus looks 
promising. Unfortunately, however, the photometric variability of late-type 
M-dwarfs complicates detection and requires us to carry out a more realistic 
appraisal of the detectability of the planets that form in our simulations. 

Fully convective stars near the bottom of the main sequence often display 
sign ificant magnetospheric a ctivity and accompanying photometric variabil- 



ity (IRockenfeller et al. 



20061 ). This is demonstrated by the variability of M- 



dwarfs, e.g. Proxima Centauri, which exhibits X-ray flare activity with possible 
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attendant photometric variation. In general, this variabihty can be modeled as 
longer-period correlated or "red" noise. Red noise levels in these stars slows the 
noise attenuation of phase- folded photometry from 1 / V Ndata to ~ 1 / \/N transits 



(jPont et al. 



20061 ). The half-life of the red noise in this approximation is four 
times the transit periodicity (days to years) rather than four times the pho- 
tometric cadence (seconds to minutes). Due to this extremely slow error con- 
vergence, the utility of rapidly taking multiple exposures is greatly reduced. 
It is thus likely that a successful observational campaign will need to attain 
cadence and signal-to-noise levels that allow the detection of a transit to high 
confidence on the basis of a single event. The candidate transit then triggers a 



relatively intense follow-up campaign that ideally 
i s in p rogress. The MEarth survey described in 



jegins while the initial event 



Nutzman and Charbonneau 



(120081 ) provides an example of this type of detection strategy and the calcu- 



lations that follow are based on their scenario. 

Our approach requires an estimate of stellar photometric noise over a time 
scale of hours. As a concrete example, we t ake o ur noise model from photo- 



metric observations of GJ 436 flGillon et al 



20071 ) as obtained using the Wise 



Im telescope on April 24th, 2007. This data set illustrates, to a first approxi- 
mation, the sort of data that might be regularly obtained by a ground-based 
survey using dedicated small telescopes. The Gillon et al. photometry was ob- 
tained during an 2.9-hour observation centered on the transit of Gliese 436 b. 
Adopting only the photometry outside of transit yields 177 photometric points 
over 2.12 hours with a standard deviation of 0.39%. No previous treatment has 
been applied to the data aside from the standard reduction and differential 
photometry. We process this data to produce a photometric red noise signal 
that can be apphed to generate Mon t e-Car lo data as follows. 

— 1 1 — I 

First, we distribute the iGillon et al.l (l2007i ) photometric data into 10-minute 
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Gillon et al. 



(120071 ) photom- 



bins, and use the photometric scatter within these bins to cahbrate a model 
for the white noise contribution. We generate synthetic 177-point datasets of 
pure Gaussian noise and calculate the standard-deviation within bins of 14 
data points. We repeat this procedure until the restricted standard-deviation 
converges and the difference from the known noise amplitude supplies our cor- 
rection factor for the white-noise amplitude giving a^hue = 0.36%. We then use 
the Scargle periodogram of the remaining longer-period noise (the smoothed 
photometry) at periods from 10 minutes to 12 hours to generate synthetic 
177-point datasets. We sum one hundred cosines of random phases with the 
amplitudes given by the periodogram to the white-noise. We then iteratively 
adjust the amplitude of the long-period noise unt il the overall noise a mplitude 
of the synthetic datasets converges to that of the 
etry. 

Our estimate of the longer-period red noise, ared = 0.15%, is intended to be 
independent of the noise introduced by the telescope and is assumed to arise 
from intrinsic stellar variability on short timescales. Red noise will limit the 
sizes, and hence the masses of the planets to which our observational campaign 
is sensitive. A search strategy that forgoes folding must be able to detect a 
transit with sufficient signal-to-noise to minimize false-alarms and expensive 
large-telescope follow-up to an acceptable level. Thus, even with zero white- 
noise, we are limited to stars that are small enough to show a transit from one 
of our planets that rises significantly above the red noise level. For example, 
if we demand 4cr confidence in the reality of a given transit signal, we could 
survey only stars small enough to yield transit depths of at least ~ 0.6%. 
The red noise amplitude will vary with time, and from star to star. Hence the 
estimate of ared = 0.15% adopted here should be understood to provide only 
a starting point for a more detailed investigation. 
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Obviously, our strategy benefits from a ready supply of small nearby stars. 
To quantify the available census, w e adopt the Lepine-Shara Proper Motion 



Catalog - North (ILepine and Shara ^ 



pie of the proper motion catalog ( jLepine 



20051) . We b egin with the local subsam- 



20051 ) consisting of 4131 dwarfs. 



subgiants and giants locate d with in 33pc of the Sun. Following the lead of 



Nutzman and Charbonneaul (120081 ) . we use color cuts to restrict the list to 



2401 local M-dwarfs by requiring that H-K < 0.7, J- H > 0.12, J-K > 0.7 



mass-luminosity relation and then we use the 



' 1 ' 1 

using the iDelfosse et al.l ( 


2000) 


Bavless and Orosz 


( 


2006 


) em- 



pirical mass-radius relation to approximate the stellar radii. Given the stel- 
lar radii, we further restrict the sample to harbor only stars with radii small 
enough to allow detectable transit depths (above a threshold of TDcru = 4cr = 
0.6%). We thus adopt a transiting planetary radius of 1.2R^, and require this 
size planet's transit signal to be equal to or greater than TDcrit- This adopted 
planetary radius is chosen so that we will only exclude from our target lists 
those stars which are so large that detection of any of our formed planets is 
unlikely. In effect, this restricts our target list to only those late-type M-dwarfs 
of radius i?* ~ O.12i?0. This profoundly limits the sample, leaving a prime 
catalog of 169 local, very late-type M-dwarfs. These stars offer the single best 
opportunity to detect a habitable transiting planet from the ground. 

We need to monitor M dwarfs with repeated visits at a cadence that is slightly 
more frequent than the expected transit duration. We must therefore evalu- 
ate the signal-to-noise levels f or any star-telescope combination a s a function 



of integration time. Follo wing 



Nutzman and Charbonneaul (120081 ) we use the 



bolometric corrections of iLeggett et al. 



(I2OOOI ) to derive the stellar luminosi- 
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ties. Luminosities and radii yield values of the stellar effective temperatures, 
Teff. Our mass and radii estimates then provide surface gravity (log (7) esti- 
mates for our stars. With Tp/-/- and lo gq, we can specify a synthetic stellar 



spectra, /(A) (IHauschildt et al. 



19991 ). When combined with a target's dis- 



tance, d, this synthetic spectra provides an estimate of the incident stellar 
flux 



Fi\) = /(A) 



-R* \ A 



d J he 



(2) 



where F{X) has units of photons/cm^ /sec/ A. With these assumptions, the 
number of photons received, Nph, from a given target in exposure of length t 
on a telescope of diameter D, is 



Nsph = t X 7r(D/2)2 x J F(A)T(A)dA 



(3) 



where T(X) is the telescope's transmi s sion 



mimics 



unction. Our transmission function 



Nutzman and Charbonneaul (120081 ) 'i-fz' filter that opens at around 



700nm and stays open for longer wavelengths, combined with the quantum effi- 
ciency of a typical commercially available CCD (which falls redward of 800nm 
and completely closes off by lOOOnm). We enforce an additional throughput 
loss of 50% to account for other unaddressed potential effects. 
With a photon count, we can obtain an accurate estimate of the photometric 
signal produced by a given target. We approximate the signal-to-noise ratio 
of an exposure as: 



SNR 



ph 



Nph + N'^hi.^lcint + ^led) + npi:,{Nskytobs + Ndarktobs + N'^ead. 



= (4) 



where ared = 0.15% as sta t ed ab ove, and agcmt is given by the scintillation 



expression of iDravins et al 



(1l998l ). The quantity npi^ represents the expected 
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number of pixels covered by the stellar image. Ngky — 10~^ e~pixel~^sec~^ 
is an estimate for the sky brightness. N^ark = 1-4 • 10~^ e~pixel~-'^sec~^ is an 
estimate of the dark-current noise based on the Nickel 1-meter telescope's di- 
rect imaging camera CCD at Mount Hamilton. Nread — 11.7e~pixel~^ is the 
amplifier-generated readout noise, again based on the Nickel 1-metcr's CCD 
specifications. For reasonable ground-based exposure times the readout noise 
dominates over the dark-current and sky, contributing ~ 0.4 mmag noise. 
Thus, for red noise levels below 0.4 mmag the readout noise of our camera 
will begin to dominate the SNR value. For fixed telescope size and exposure 
time, the required SNR yields a tradeoff between spurious transit-like noise 
and missed transits. Our goal is to optimize the choices of telescope diameter, 
exposure time, and SNR value for recovery of our simulated planetary popu- 
lation. Our exposures must be a few times shorter than the expected transit 
duration. For close-in terrestrial planets whose transit durations are less than 
60 minutes, this limits the exposure time to less than about 20 minutes. Given 
this maximum exposure time, we can estimate iV ~ 21 observations per 7-hour 
night. Three targets on a given night yields ~ 7 photometric observations per 
star. We thus also require that noise be discernible from a transit signal at the 
4.3(7 level, allowing ~ 99.8% confidence in any given (~ 7 point) lightcurve. 
We adopt telescope sizes of 1 meter in keeping a realistic expectable financial 
investment for a dedicated array. 

Figure 7 shows the number of stars that are bright enough to allow observation 
within 20-30 minute exposures for varying stringencies of the required SNR. 
Note that there are few stars which allow a statistically significant detection. 
For example, if we adopt a detection criterion of 3cr, the number of avail- 
able target stars is unacceptably low (~ 10). Thus, for correlated noise levels 
c^red ~ 0.15%, the prospect of detecting transiting Earth-sized planets with 
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an array of meter-class ground based telescopes is challenging. Our Monte- 
Carlo detection simulations therefore examine the effects of different levels of 
correlated noise. 



Our observability simulations extend the work of ISeagroves et al.l (120031 ) and 
specific details of the computational procedure can be found in that work. A 
Monte Carlo simulation is initialized with the planetary populations (drawn 
from our simulations), a list of observers, and a stellar target list (the low-mass 
stars defined by TD > TDcru)- Each observer has an associated location (lat- 
itude, longitude and altitude) and weather (average fraction of clear/cloudy 
nights per year). Each target has an associated position (R.A. and Dec), stel- 
lar mass, and a photon flux over the wavelength region of interest. The Monte 
Carlo routine starts by assigning which targets will host transiting planets in 
the simulation. The program assigns planetary systems (the end-state of our 
N-body simulations) to a fraction of these target stars and then orients those 
planetary systems randomly with respect to the line-of-sight to Earth. As- 
signed planets maintain their orbital periods, and hence assume a semi-major 
axis in concordance with the parent star mass and Kepler's third law. The sim- 
ulation proceeds through an observing campaign night by night. Every night, 
the simulation must determine where the weather is favorable: while season- 

:he fraction of 
The night's 

weather for each location is determined by drawing from these probabilities. 
If clear or cloudy, all observers common to the location will be affected identi- 
cally. If the weather is determined to be partly cloudy in a given location, then 
it is possible that some observers at that location will be able to observe while 



based weather patterns have not been figured into the model, t 
clear and cloudy nights at each observing location is knowij^ 



^ http:/ /www. ncdc.noaa.gov/oa/climate/online/ccd/cldy.html 
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others will not. For each observer at a location with favorable weather, the Ju- 
lian dates (JD) of sunrise and sunset at the observer's location are calculated 
for the current date in the campaign. The JD of sunset, the position of each 
target, and the observer's location are used to calculate air masses for each 
target in the target list. Any targets that pass the airmass cutoff (sec(z) = 2.5 
in our simulations) at sunset are then checked to ensure they will pass the 
airmass limit for at least 2 hr. Each observer will then construct a target-list 
from those targets that have passed the airmass cutoffs and whose required 
exposure times are short enough to allow multiple targets to be observed every 
hour. The observer then cycles through this observing list every hour until one 
of the targets fails an airmass test, at which point it is dropped and a new 
target is added to the observing list. On each visit, an observer integrates long 
enough to reach a 5a signal-to- noise ratio before slewing to the next target (1 
minute of combined slew and readout time is added between each target in 
the observing schedule). This level of signal-to-noise is higher than that used 
to signal a detection so that we are able to detect transits more quickly and 
reliably. Once a detection has been made, we anticipate extensive follow up 
observations by larger telescopes, however in this preliminary study we do not 
modify our observing plan as a result of transit detections. The list is repeated 
every hour so as to catch the majority of transits that might occur over the 
course of the night. In the event that one of the targets hosts planets, and 
that one of the planets happens to transit during the night's observations, 
the line-of-sight separation between the planet and its host star is calculated. 
Using this separation, as well as the radii of the transiting planet and its host 
star, the photometri c trans it depth is determined using the approximate for- 



mulae of 



Ohta et al. 



(120051 ). This process is repeated for each observer in the 



observer list, and the aggregate process constitutes 1 night of the campaign. 
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The simulation then increments its internal calendar by 1 day and repeats the 
procedure for the overall duration of the campaign. Using this rotating obser- 
vation method rather than continual monitoring of individual targets allows 
us to monitor several stars simultaneously with only minimal loss of coverage: 
the typical transit length is one hour, so one observation every hour should 
catch any transits that occur. In comparison, the typical orbital period is be- 
tween 10 and 60 days, so if we wished to monitor individual stars continually 
we should do it at least over this time range. This severely limits the rate at 
which we could study our targets, extending the necessary time to discovery 
by many years. 



The results of the Monte-Carlo simulation depend on the number and distri- 
bution of observers, the stringency of the detection criterion (SNR- value), and 
the adopted level of correlated noise. We set the fraction of stars with planets 
to 100%, this result is of course downwardly scalable. As an illustrative case 
of a viable search strategy, we distribute 10 observers in longitude around 
the Northern hemisphere. We choose this distribution of observer locations 
envisionin g the utilization of a p re-existing infrastructure of small telescope 



observers (ISeagroves et al. 



20031 ). and we note that such a distribution could 



be altered more observers at fewer sites without a significant impact on the ca- 
dence or quality of the results. We choose this distribution of observer locations 
envisionin g the utilization of a p re-existing infrastructure of small telescope 



observers (ISeagroves et al 



20031 ). and we note that such a distribution could 



be altered more observers at fewer sites without a significant impact on the 
cadence or quality of the results. Each observer operates a 1-meter telescope 
and signals a potential transit whenever a photometric measurement falls A.Sa 
below the night's average. This stringent detection criterion of 4.3cr was cho- 
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sen to reduce false-positive occurrence to acceptable levels. Due to the highly 
non-Gaussian nature of the false-positive rate, we iterated to determine this 
optimal value of A.Sa where we were able to reduce the false-positive rate to 
(on average) 19 over the entire two-year campaign while still allowing rea- 
sonable integration times. We vary the amplitude of the correlated red noise 
from 1.5 mmag to zero, and for each red noise level we repeat the model ob- 
serving procedure several (5-10) times. The averaged results stemming from 
135 Monte-Carlo simulations are presented in Fig. 8. A ground-based transit 
search of the above specifications is able to detect Earth-sized planets with 
acceptable levels of false-positives when red noise is controlled to the level of 
0.45 mmag (30% of the adopted value). Over the course of a two year cam- 
paign the observers produce approximately 9000 photometric light-curves. Of 
these light-curves, on average 5.5 contained true transits and 1.3 of these were 
detected at the reduced level of correlated noise (0.45 mmag). At the end of 
these campaigns there are another 1.5 transiting planets that could be detected 
but haven't yet, thus a campaign that utilized more observers or an extended 
search would be expected to detect a total of ~ 2.8 transiting planets given 
the reduced level of red noise. If we utilize the h = scaling occurrence frac- 
tion of / = 28%. We obtain an estimate of 0.8 detections. Our results are 



tabulated based on planetary radii drawn from the results of 



Valencia et al. 



(120071 ) that assume ice-mass fractions of 0%, 50%, and 100%. This allows us to 
cover a reasonably expected range of possible planet compositions and sizes, 
the resulting planetary radii, number of planets detected under said observing 
campaign and the number of potentially detectable planets under this observ- 
ing configuration are presented in Table 2. These other simulated observing 
campaigns are only different in the planetary radii used, all other parameters 
were held constant to facilitate comparison. We find that the number of de- 
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tections is small regardless of the ice-mass fraction of the planets, similarly 
the number of potential detections is uniformly small, making a statistically 
significant null-result elusive. Irrespective of the red noise level, given our SNR 
detection threshold of 4.3o", there were an average of 19 spurious false-positive 
'transits' reported over each of the 2-year observing campaigns. 



Transit detection from space has been a topic of discussi on following the direc t 



detection of HD 209458b by the Spizter Space Telescope (IDeming et al 



as we ll as the constraint on its albedo by the MOST telescope (iRowe et al. 



2005) 



20061 ). To examine the potential detection of our simulated planets using a 
dedicated space-based mission, we ran a simulation meant to emulate the 
general characteristics of the proposed TESS satelhte. As before, we use our 
suite of planet-forming simulations to populate the target stars. The TESS 
satellite is designed to have six small (~ 13cm) wide-field cameras observing 
stars in the direction of the anti-solar point, with a total field of view of ~ 2000 
square degrees. This setup will yield ~ 36 days of continuous photometry for 
each star over it's two year survey. Using a simplified version of our Monte- 
Carlo code, with a fixed exposure time of 1 minutes, we simulated observations 
of the 1385 best candidate stars in the LSPM-North catalog. This sample of 
stars was determined as described above but with the stellar size cut {TDcru) 
based on the reduced noise level due to the lack of atmospheric scintillation. 
We then analyzed the resulting lightcurves by sliding a characteristic transit 
template along the lightcurve and for each time- alignment calculating a 
value between the lightcurve and the template. We then divide this value 
by the of ^ iiull hypothesis (no transit). We fold this figure of merit time- 
sequence along the possible planetary periods, and search for periodic signals. 
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Of the 144 stars that had transits, 31 were clearly detectable while the rest are 
too small relative to the stellar (red) noise to be detected. If we assume that 
the Southern sky would allow another 31 detections and use an occurrence 
fraction of 28%, we find that such a mission would allow the detection of ~ 17 
of the planets in our simulations after the two year survey. 



5 Summary 



It is generally agreed that if terrestrial planet formation around low-mass M- 
dwarfs is a scaled down version of solar system formation, then planets in 



the habitable zones of low-mass red dwarfs will be sma^ 
inhospitable (see, e.g. 



Raymond et al. 



2007 



Lissauer 



1, dry , and thoroughly 



20071 ). Nature, on the 



other hand, tends to defy expectations. Our tack here has been to investi- 
gate an alternative, optimistic scenario for terrestrial planet formation around 
nearby low mass stars, which assumes that disk-mass is largely independent of 
stellar mass. In this scenario, the planets that emerge are more akin to the Jo- 
vian satellites rather than Mercury, Venus, Earth and Mars. Our motivation 
for this investigation stems from the fact that our hypothesis is potentially 
testable, and at a relatively low cost. 

Our approach couples end-to-end simulations of both the in-situ formation 
and the transit-detection of terrestrial mass planets around local, late-type 
M-dwarf stars. Our formation simulations followed the accretional growth of 
planetary embryos, around an = 0.12 Mq star, from the onset of oligarchic 
growth {Memb ~ 3 X 10^^ g) using a semi-analytic method, through the stochas- 
tic late-phases of growth using a full 3-D Nbody simulation. These simulations 
resulted in systems of several (3-5) planets with masses comparable to Earth in 



30 



stable orbits in or near the star's habitable zone as well as 1-2 planets per sys- 
tem with masses comparable to Mars. We found that these results were largely 
insensitive to the presence or absence of an exterior perturbing Neptune-mass 
planet, but detectable planets are more likely to occur if an external perturber 

is present. 

Photometric red noise makes transit detection of these small planets challeng- 
ing. We have explored the impact of varying levels of red noise on the potential 
utility of dedicated ground-based telescope arrays. Our method consisted of 
a detailed Monte-Carlo transit search code implementing realistic weather, 
stellar flux, telescope and M-dwarf noise models. We find that a formation 
frequency of ~ 28% and nightly stellar noise levels in the vicinity of cr ~ 0.45 
mmag allow a longitudinally distributed array of ten dedicated meter-class 
telescopes carrying out a targeted search to detect, on average, 0.3 transiting 
terrestrial planets within two years of operation and another 0.5 potentially 
detectable planets. However, a sateUite based survey similar to TESS is capa- 
ble of detecting ~ 17 planets after it's two year survey. 

Time and again over the past decade, nature has surprised planet hunters with 
the sheer diversity of planetary systems. Experience has shown that planet for- 
mation under a variety of conditions can be more efficient than the example 
of our own system would suggest. Our approach in this paper has been to 
investigate a planet building model that lies at the optimistic, yet still justifi- 
able range of formation scenarios. Should our picture prove correct, then the 
prospects for discovery of truly earthlike, alarmingly nearby, and potentially 
habitable planets may fie close at hand. 
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Final Planetary Configuration Statistics 













Simulation Group 


N Planet 


Mpianet 


Mass Lost 




A No Ice-giant 


4.4 


0.79 


2% 


1.4 


B Ice-giant present 


3.8 


0.68 


18% 


0.2 



Table 1 

"Planet" indicates masses above 0.05M®. Npianet is the average number of planets 
per system. Mpi^net is the median mass of planets per system. Np is the average 
number of planetary embryos that survived per system. 

Comparison of Planetary Compositions, Radii and Detectability 



Ice 


Rp 


^detect 


■^possible 


100% 


1.46 


0.4 


0.6 


50% 


1.20 


0.3 


0.5 


0% 


0.93 


0.1 


0.3 


50%* 


0.68 


0.1 


0.1 


50%** 


0.24 


0.0 


0.0 



Table 2 

Columns indicate the planetary ice- mass-fraction (Valencia et al., 2007), median 
resulting planetary radius (in Earth radii), the number of planets detected in our 
observing scenario after the two year simulated observing campaign, and the number 
of additional planets potentially detectable given the levels of rednoise present, 
respectively. Assuming an occurrence fraction of 28%. 

The bottom two rows respresent our low disk-mass simulations with normalizations 
of * 7g/cm^ and ** 0.84 g/cm^. 
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Fig. 1. Beginning/end (thin/thick lines) surface density profiles of the planetary 
embryos and planetesimals (solid and dashed lines, respectively) from the semi-an- 
alytic model. The semi-analytic model runs until ~ 150 years. 
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Fig. 2. Evolution of Simulation A-01 (no ice-giant present) ending with planets of 
mass 0.86, 0.79, 1.03, 0.36 and 0.22 (in order of increasing semi-major axis) 
after 10^ years. The size of the symbol is proportional to the radius of the body. 
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Fig. 3. Evolution of Simulation B-01 (ice-giant at 0.6AU) ending with planets of 
mass 0.60, 0.75, 0.83, and 0.41 (in order of increasing semi-major axis) after 
10^ years. The size of the symbol is proportional to the radius of the body. 
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Fig. 4. Final planet masses using a non-uniform initial mass distribution. Group 
A is represented by the solid line (no ice-giant present). Group B is represented by 
the dotted line (ice-giant at 0.6AU). 
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Fig. 5. The amount of mass located in each semi-major axis bin. Group A is 

represented by the sohd hne (no ice-giant present). Group B is represented by the 
dotted hne (ice-giant at 0.6AU). 
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Fig. 6. Probabihty of detecting a transit of a planet around its host star as a 
function of the telescope's sensitivity. Dotted hne shows the simulations with an 
ice-giant, solid hne shows the simulations without. An approximate sensitivity of 
sub-meter class telescopes is shown for comparison (~ 3 milli- magnitudes) . 
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Fig. 7. Number of stars which could yield detectable transit signals in a 20 (30) 
minute exposure as a function of signal-to-noise requirement, shown by the solid 
(dashed) line. 
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Fig. 8. Number of transit detection after two years of observation by a dedicated 
array of meter-class telescopes, as a function of red noise amplitude. Scaled to an 
occurence fraction of 28%. 
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